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Abstract. We discuss the existence of a dividing shell separating expanding and collapsing regions in spherically symmetric 
solutions with pressure. We obtain gauge invariant conditions relating not only the intrinsic spatial curvature of the shells to 
the ADM mass, but also a function of the pressure which we introduce that generalises the Tolman-Oppenheimer-Volkoff 
equilibrium condition, in the framework of a 3+1 spacetime splitting. We consider the particular case of a Lemaitre-Tolman- 
Bondi dust models with a cosmological constant (a A-CDM model) as an example of our results. 
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INTRODUCTION 

Cosmological structure formation assumes the collapse of inhomogeneities, via gravitational instability, into "bound" 
structures, with the underlying idea that they depart from the cosmological expansion. This approach usually models 
overdensities within closed patches embedded in Friedman backgrounds, in particular in the spherical collapse 
included in the Press & Schechter scheme [1]. Birkhoff's theorem is often invoked to claim that the evolution of 
the overdensities is independent |0], while, rigourously, it only applies to asymptotically fiat spacetimes from the 
cosmic expansion [I3[] . 

In this work, we define general conditions for the existence of a shell separating collapse from expansion, and 
will illustrate our results with a simple example of inhomogeneous A-CDM models. Our work differs from previous 
approaches (see e.g. |4j,|5[]) since it does not involve spacetime matchings or metric perturbations. Instead, we adopt the 
Generalised Painleve-Gullstrand (hereafter GPG) formalism used in Lasky & Lun [6], which involves a 3 + 1 splitting 
(ADM) and the consideration of gauge invariants kinematic quantities 071]. We then define general conditions for the 
existence of a shell separating contraction from expansion, before proposing particular examples of inhomogeneous 
A-CDM models implemented as Lemaitre-Tolman-Bondi models with dust and a cosmological constant. Finally, we 
perform a dynamical study of such models. 



ADM APPROACH TO LTB MODELS IN GPG SYSTEM 

We consider a spherically symmetric Generalised Lemaitre-Tolman-Bondi metric to include pressure. Performing an 
ADM 3+1 splitting in the GPG coordinates |0] , the metric reads 

ds 2 = -a(t,r) 2 dt 2 + 1+ * , (P(t,r)dt + dr) 2 + r 2 d£l 2 , (1) 

where a is a lapse function, j3 is a shift, and E > — 1 is a curvature-energy characterising the curvature of the spatial 
surfaces orthogonal to the direction n a = (—a, 0,0,0) of the flow. For a perfect fluid, the projected Bianchi identities 
r, a fl = yield the energy density conservation equation after projecting along the flow n b , and the Euler equation 
projecting orthogonally to n b . Using the projection h%, and the Lie derivative of the density p along the flow Jzf„p, we 
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FIGURE 1. Kinematic analysis of motion in the pseudo-potential V. Depending on E relative to £/,„,, for a given shell of constant 
M and E, the fate of the shell is either to remain bound (£< < En m ) or to escape and cosmologically expand (£> > Eu m ). There 
exists a critical behaviour where the shell will forever expand, but within a finite, bound radius (E = £;, m , r < r/ im ). 



have: 

n b T§. a = -J? n p-(p+P)© = 0, hXc =0 => P' = - (P + P) ^ , (2) 

where P is the pressure, the radial derivatives are denoted by a prime, ', the time derivatives are represented by a dot, ', 
and the expansion is ©. 

Introducing the ADM (also called Misner-Sharp) mass 

M = r 2 {l+E)(\na)' -4-nPr 3 + ^Ar 3 + r 2 ^„ , (3) 

where A > is the cosmological constant, we write Einstein's field equations (EFEs) as Lie derivatives along the flow 
and the radial evolution 

X,M =AnPr 2 V ^M=j3(M'+ AnPr 2 ) , (5) 

J3\ 2 M 1 , 

tL\ =E + 2- + -Ar 2 . (6) 
a J r 3 

The system becomes then closed when an equation of state is supplied. The A term can be absorbed in M, E and P, 
and many fluids mass equations written with each component using the ^ term for the overall sum of the masses. 

SEPARATING COLLAPSE FROM EXPANSION 

To characterise the separation of expansion from collapse, we shall use GLTB coordinates in which the flow direction 
reduces to dj by choosing j3 = — r. We thus have 

ds 2 = -a(T,R) 2 (d T t) 2 dT 2 + ^* r ) dR 2 + r 2 dQ 2 (7) 
^ ; v ' l+E{T,R) 

i-i „ 1 +E , ( r \ 2 Ml, 



M =P4nPr 2 ,Er = 2j3 P', =£ + 2— + -Ar 2 . (8) 

where t and r are functions of the coordinates T and R. There are two situations one should consider in parallel. 
On the one hand, we look for the gauge invariant expansion ©, defined as n a - a , since our goal is to separate an inner 
collapsing, spherical region from the outer expanding universe. On the other hand, the total ADM mass of the spherical 




FIGURE 2. Plots of log(-£,, m ) - log(fi) and log(-E) - log(tf) for a) NFW H with background and an E profile given by 
E(R) = —4E m j„ [J^j [1 — ^> for E m [ n = — 1 + e~ 10 and r\ = e 9 , and for b)power law density without cusp + background and 
Hubble initial flow. 



region that departs from the expansion flow should be conserved. This is indeed suggested by the dust case where that 
happens for every shell. 

Denoting with * an evaluation at the dividing shell, we find the following relation between the expansion © and the 
shear a, in that at the * shell 

= 0^ JSVI* =0 when Jz?„M(f ,r*(f)) = 0. (9) 

Moreover the generalized Friedman constraint 

(3) fl + ^© 2 =6a 2 + l6np + 2A (10) 

tells us that the vanishing of © only happens in regions of positive 3-curvature ^R. On the other hand, if we demand 
that the separating shell has a dust-like vanishing mass/energy flow, i.e., has a conserved ADM mass along n": 

Vf, Jz?„M(f,rJf)) =0 +>Vt,E = -2 — --Arl < 0. (11) 

r* 3 

We remark that M refers to the total ADM mass, thus including a cosmological constant. 

Taking into account that the equilibrium of static spherical configurations requires the satisfaction of the Tolman- 
Oppenheimer-Volkoff equation of state [9], we are led to define a generalized gTOV function 

which reduces to the usual TOV equation when it vanishes on the ★-shell. 

The radial behaviour of the *-shell is then similar to a turnaround shell. Indeed = — W*- leads to a null radial 
velocity r* = while its acceleration reveals the importance of the gTOV parameter 

?GLTB,*=-cc 2 gTOV ic . (13) 

The local staticity of the ★-shell is then shown to be equivalent to having a local TOV equation on this limit shell 

gTOV^ = ^%}r* = _Sf„ (0 + 3aY =0 <^ — P' = 

p+P 
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FIGURE 3. Closed and flat backgrounds with outwards escaping shells: protected inner shells. 



AN EXAMPLE: ACDM 



A simple illustration of our result is given by the case of dust with a A. There is then no pressure gradients and M^ m 
and E are conserved, i.e. a = 1, which simplifies the analysis and allow us to perform a kinematic study (see Fig. 
Q} per shell of Eq. |8}c). An effective potential is defined by V(r) = — ™ — jr 2 , so that when r = we have indeed 

V(r)=E: 



M 1 



-Ar z +E, 



with f - 



M 



(15) 



(3M)3 A3, which only depends on A and M(R). 



For each shell there is a virtual static state at r/ lm — , — 

We also recall that, in this model, gTOV = — = —r. 

We then need to choose initial conditions p,-, which set the profile, and v,, which set the £, profile. The inter- 
section of En m with Ei will be static at turnaround, hence defining our limit shell. We use two sets of cosmologically 
motivated initial conditions: one with p ( - as a well known "universal" CDM halo density profile from simulations, 
refered to as NFW profile [8] and as a parabola (Fig.[2]-a), and the second, with a cuspless power law as p,- and a 
Hubble flow for V; (Fig. 0-b). Explorations of all general cases for initial conditions yields a split between the sep- 
aration of inner shells from outer influence, which is guaranteed for ever expanding backgrounds (thus unless the 
background recollapses, see Figs. [3] and |4j and that of outer shells from inner influence, which occur for open and 
some fiat backgrounds (see Figs. |4j». 



CONCLUSIONS 

Using non-singular, generalized Painleve-Gullstrand coordinate formulation of the ADM spherically symmetric, 
perfect fluid system [6] we have shown ifloifTTll that the existence of shells locally separating between inner collapsing 
and outer expanding regions, is governed by the condition that the combination of expansion scalar and shear 6 + 3a 
should vanish on the shell. The ADM mass of the shell is then conserved. This condition requires that the separating 
shell must be located in an elliptic {E < 0) region. Moreover, for that shell to exist over time, we have shown 
that the TOV equation must be locally satisfied. We argue in some cases that this local condition is global in a 
cosmological context (FLRW match at radial asymptote). Given appropriate initial conditions, this translates into 
global separations between an expanding outer region and an eventually collapsing inner region. We present simple but 
physically interesting illustrations of the results, a model of Lemaitre-Tolman-Bondi dust with A representing spherical 
perturbations in a ACDM model with two different initial sets of cosmologically interesting conditions consistent with 
known phenomenological constraints [fTol [Til and Refs. therein]: an NFW density profile with a simple curvature 




FIGURE 4. Flat, inner shells free, and open backgrounds: two limit shells 



profile going from bound to unbound conditions, and a non cuspy power law fluctuation with initial Hubble flow. We 
show, for these models, the existence of a global separation. We argue that these shells are trapped matter surfaces 
ifToL [Till and that, therefore they separate domains of influence of cosmic expansion from local domains of matter 
dynamics and can hold the place, in structure formation studies, of the incorrect invocation of Birkhoff 's theorem. 
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